NORDITA-95/3 P 
January 1994 
hep-th / yymmddd 



THE GROSS-NEVEU MODEL ON THE 

LIGHT CONE 

Igor Pesando []|] 

NORDITA 
Blegdamsvej 17, DK-2100 Copenhagen 
Denmark 



Abstract 

We consider the (massive) Gross-Neveu model using the hght cone 
quantization where we solve the constraints explicitly. We show that 
the vacuum is trivial and that the quantization fails when m = 0. We 
discuss how the running coupling constant emerges as a pure normal 
ordering effect in this context and the bound state equation. 
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1 Introduction. 



One of the major outstanding problems of the contemporary physics is to 
find a way to compute observables in strongly interacting field theories, of 
which the theory for the strong interactions QCD4 is the best examples. 
Many methods had been used but for now the most promising are lattice 
calculations and light cone field theory: the latter has received a lot of at- 
tention during the last period. The advantages of the light cone are that the 
maximum number of Poincare generators become kinematical (at = 0) 
and that it is very easy to write a formula for the mass spectrum. Until now 
the approach has been successfully used to solve a variety of 2D problems 
such as large N QCD2 with vector matter (0) and adjoint matter ([§]), and 
applied to (0^)1+1 (i) and QED^ ([|). 

A different approach based on the renormalization group has been advo- 
cated and is developing by Wilson and collaborators (see for example (|^) 

In view of a better understanding of the light cone approach we consider 
the first non trivial, i.e. non superrenormalizable model in 2D, the Gross- 
Neveu model. The (massive) Gross-Neveu model was previously treated in 
the infinite momentum frame (see for example (0)) and in the light cone 
without solving explicitly the constraint and assuming the triviality of the 
vacuum (0). Instead we will solve the constraints explicitly and we will prove 
that the vacuum is trivial but this approach fails when m = 0, nevertheless 
it yields the running coupling constant as a pure normal ordering effect and 
the range of stability of the theory in a very simple way. 

2 Gross-Neveu model in the light cone. 

The lagrangian of the (massive) Gross-Neveu model ([jl[) is given by 

C = ■ {i ^ -m)ij + j^ii; ■ ifjf (2.1) 

that can be explicitly written in the light cone as[] 
^ Conventions. 

^ ^ ^{x° ± x^) A''B^ = AoBo-AiBi=A+B^+A+B^ e"^ = -e+- = 1 
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(2.2) 

where = [ip^) = if)* with i = 1 . . . N . As it is usual in the hght cone 
approach primary constraints are given by the classical equation of motion 
for the nonpropagating fields 'i/'L 

- iV^d-i^l - mi)\ + 2|^(V'+ ■ ^A- + ^- ■ 4^+) = (2-3) 



and -i/^! 



iV2d^^'_ - rm/j^ + 2^(^+ ■ i/j- + ifj- ■ ^+ = (2.4) 
Using these constraints we can rewrite the lagrangian ( |2.2| ) as 

£' = tV2i;+ ■ - ^(^_ . + V;+ ■ ^_) (2.5) 

where -0- is to be seen as a functional of ip^. From the previous effective 
lagrangian we get the translation generators as 

TTi f — — 

P~ = Y y dx~ ■?/'+■ -0_ + ■0- ■ 

= 1^2 j dx-ip+ ■ d_il)+ (2.6) 

These generators are hermitian because we started from a real lagrangian, 

i.e. with the explicit {ip- ^ ■0); would we not have used such a real lagrangian 
we would not have got hermitian generators. 

/0\/2\ /00\ /^01\ /^Ol 

^+ = (, j = (, V2 j = (, 1 j ^^ = (,-10 

/ 1 \ p 1±75 
75 = -7o7i = n 1 Pr.l = 



0-1 "-"^ 2 



= d 



j2„ 
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Now we quantize with the standard Dirac brackets 

{V';(x),^U2/)}lx.+=.+ = ^5^^6{x~-y-) (2.7) 

in the hght cone box x~ G [— L, L] where we impose the standard antiperiodic 
boundary condition 

i)\{x- + 2L) = -ijiix-) ipXix- + 2L) = -^V(x-) (2.8) 

When we expand the operator -0+ in Schrodinger picture in Fourier modes 



v2 rjj 1 ^ ylL 
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eq. ( |3.4| ) is equivalent to impose the anticommutation relations 

{i^l^^lji} = 5rs5'^ (2.10) 

With in mind the idea of using a variational approach similar to that used 
by Coleman (0) to show that the perturbative vacuum is the true vacuum 
(for a certain range of parameters), we introduce the normal order A^^[. . .] 
defined by 

_ J if r G C creation operator y _Jif r E A = C annihilation operator 

^ 1 if r G ^ annihilation operator ^ 1 if rGC = ^ creation operator 

(2.11) 

where AU C = 'Si + \, Af\C = Since we require that the vacuum \A > 
be C-invariant f\ we have to impose r E A <^==^ — r G C. The choice of the 
set A is equivalent to consider as vacuum the state 

TV 

\A>^ nnv':io> (2.12) 

sG.4 i=l 



^ We define the charge conjugation as C'4'+{x^)C^^ — {[^^{x^) -<=4> CiprC^^ = 'tp-r- 
Note that the vacuum C-invariance is a consequence of the request of having a SO{2N) 
invariant vacuum. This symmetry is evident when the lagrangian is written in terms of 
Majorana fermions. 
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where |0 > is the usual free vacuum, defined as ■0-r|O >= i/jr\0 >= for 
r > 0. 

After this introductory stuff we can try to solve the constraints explicitly 
and then to write down the explicit form of the translation generators ( p.6| ) . 
Actually, we are not interested in solving eq.s ( ^]4| ) and ( |2]^ ) but our aim 
is to express the operator T = ip^ in function of the propagating 

fields ip^ and '0+. 

The formal solution of the constraint ( |2.4| ) reads 

= --^d-J^X + ^az^[^^UV^+ ■ + ^_ ■ ^+)] (2.13) 

and we take this expression to be the quantum constraint, from which we 
derive the quantum constraint for the operator ip- by hermitian conjugation. 
From this equation we derive the equation for T 

T = -^(^+ ■ azV+ - az^^+ ■ ^+) 

+ -^1^+ ■ dz'ii^+T) - dz\Ti,^) ■ (2.14) 

The solution of this equation can be given as an expansion in power of the 
bilinear ■0+(x) ■ ip+iu)] the explicit solution in the leading order in is given 
in the appendix, here instead we want to make some comments: 

• Since we quantize in a finite box we want to impose antiperiodic bound- 
ary conditions on ip^ and on ifj_ (obvious if we think to the Dirac 
procedure). 

This means that the solution ( |2.13| ) has to "propagate" the antiperi- 
odic boundary condition from ip^ to ip^ and this in turn requires that 
we define properly dz^, i-e. in such a way that ip- satisfies antiperiodic 
boundary conditions. Explicitly we have 

d-\x)=y- (2.15) 



Notice that there is a possible source of ambiguity: this is due to the 
order of ip'^ w.r.t. [ip^ip- + ip^ip^) in the second term of ( |2.13| ). 
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This ambiguity does not show up to the leading order as the exphcit 
computation reveals. 



• What happens to the solution in the limit m — 0? As it is clear 
from the eq. ( |2.13| ) the perturbative solution in g'^ vanishes. This is 
not the right answer because ip^ has to be antiperiodic. Moreover the 
hamiltonian P~ also vanishes in this limit, this means that we cannot 
quantize the theory on the light cone when m = 0. If we still insist in 
finding a solution for the ■0_ at m = we have to introduce further 
constraints on the theory and as the = 1 case shows, we have to 
quantize the values of the charge J^idx~ 

As we are interested in the mass spectrum in the large N limit, we need 
not compute the exact solution of ( |2.14| ) but we can extract the leading 
contributions in j^. More precisely we want to compute the mass spectrum 
of the particles and of the "mesonic states" hence we need compute P~ in 
the approximation 

P- = Ni2L)Lo+ [ Li(a:;xi,yi)iV^[^+(xi)-^+(yi)] 

Jx,xi,yi 

iV Jx,Xi,yi 

+higher orders in i/jip 
= P^Q^ + P(l^ + P(^2) + higher orders in i/jip (2-16) 

where all the L are of 0(1). The explicit computation yields: 

P- = N(2Lf-f^' MIL^^MA^^.^,] 

^ ' 4 g^T.c{Q) - 1 2tx ^ r ^ ^ 

2^ dri+r2,si+S2 9 + Si) ■ 



jn 11 ri,r2,si,S2 



(-- 




(-- 




Vri 




\r2 





^r2 rS2 



along with 



(2.17) 



= -^7 E ^ - 7 E ^^Ai^r ■ A] (2.18) 
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where we defined the physical mass of the particles i/jrlA > (r G C)0 as 



m 



2 



and the running coupling constant 



^ = r^^TTTT^ 7T? (2-19) 



1 - g^T.c{n 

where 

2n r T + n 



= ^1 IWt (2-20) 



Sc(n) = iE^ (2.21) 



and 



S9nc{r) = \^_\ \W (2.22) 



We notice that when we normal order w.r.t. the free vacuum the running 
coupling constant can be written in the thermodynamic limit, i.e. when 
L ^ with p = ^ = const and j9 << ^ = A asQ 

9'iP) = 1 IWt = , (2-23) 

l-^2So(p) log^^5 



77 



which shows the asymptotic freedom and coincides with the usual result ([|l[). 

Before discussing the range of stability of the theory coupling constant, 
we want to point out that the vacuum expectation value of the fermionic 
condensate ■0 ■ ■0(x) is proportional to the vacuum energy density because of 
the first of {^): 

< ^ ■ i^i^) >= ^^%7^ = -iVMSc(O) (2.24) 



■^This follows from the mass formula = 2A^^[P+]iV^[P ] at the leading order. 
■* We make use of 



^ - = ln + - 



r TT A 

2 

where is the Euler constant and A is the UV cutoff. 
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which coincides with the formula (112) of ([0) in the substance. But we are 
deahng with a composite operator and this needs an independent renormal- 
ization as eq. ( ^.241 ) shows clearly since it is divergent. If we are naive and 



in analogy with the renormalization of (j) operator in a massive scalar the- 
ory in the usual lagrangian perturbation approach we take Z^^ = ^ ^ ^, 



>V dM 

which is valid in the zero momentum renormalization scheme, and we get 
< i/j ■ t/j^x) >R= —N-^^ but this is wrong because the integrated operator 
/ dx~ i/j ■ i/j is (proportional to) the translation operator P~, which has been 
renormalized such that < >= 0, hence the right renormalization is 

< ^ ■ ilj{x) >n= (2.25) 

This agrees with the zero momentum subtraction of the auxiliaty field cr (||l|), 
which is possible in the massive theory. Explicitely the effective potential 
reads 

^eff = + ^ / log(^' - ("^ - 2^^)') (2-26) 

which can be renormalized with the condition V^g'(o" = 0) = 2 and since this 
means that | =< ip ■ ip >= we recover eq. ( p.25| ). 



3 The vacuum of the theory. 

Now we want to prove that the vacuum is trivial in a certain range of and 
that outside this range the vacuum is not among the test states \A>. 
The simplest way of seeing this is to consider 

< A\No[P-]\A >= F(o)(^) - F(-)(0) (3.1) 

This quantity has to be greater than zero because the perturbative vacuum 
|0 > be the true vacuum. 

From the plot of the function P^){X) oc ("^ ~ ^(0)) from the 

bound |Sc(0)| < So(0) = ^ J2i ^ = ^max we see that two situations can arise. 
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Unstable case 



When 



9^ 



E 

2 



TT 



> 



(3.2) 



the minimum of the vacuum energy is achieved for the perturbative vacuum 
while when eq. ( p.2| ) is not satisfied, the free vacuum |0 > is not the mini- 
mum energy eigenstate of the theory, i.e. it is NOT the true vacuum. Now 
we can ask whether one of the \A > states is the true vacuum and whether 
the theory is meaningful since its ground state seems to have unbounded 
negative energy. 

In order to answer to this question and to confirm the previous coupling 
constant range we look at the mesonic spectrum and we check for tachyonic 
states. We define the following U{N) invariant operators 



1 



Mr. 



1 



ips ■ 'ipr r E C, s E A 



(3.3) 



which satisfy the commutation relations 



[M, 



risi 5 ^4s2] 




(3.4) 
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and the look for the mesonic mass spectrum. We need not the other U (N) 
invariant operators i/jr ■ ips with either r,s G ^ or r, s G C because they 
annihilate the vacuum \A > and can be always brought to act on the vacuum. 

With the help of the and M we can rewrite the leading order contri- 
bution of P~ when acting on states created by M'^ as 



NAP-] = ^ E M}J(^-U6,,.A,,, 



''^'^ ~ "'K-r.,S2-r2 (- + -](- + -] \M.,^2 (3.5) 



27r Vsi riJ \S2 r2 



It is possible to write P~ as in ( ) because we are only considering 
states generated by M^" acting on the vacuum \A >, which implies that 
we only need considering the commutation properties of P~ with the M 
operators. In particular we find 

[NA[Pn,],Ml] = ^ - - - ) Ml (3.6) 



which allows us to write F^^-, = J2s€C,r€A (s ~ r) ^l^Mrs in virtue of ( 
[31^ ). As far as P^-^ is concerned we are dropping the non leading contribu- 
tions. 

Now our task is to diagonalize the hamiltonian ( 2]5| ) 



Na[P-] = ^ E Ml,^Mr,s,,r2S2Mr2S2 (3-7) 

and to find its eigenvalues. 

To this purpose we notice that the matrix A^risi,r2S2 can be written in a 
block form where the blocks are characterized by a natural R = s — r: 



Si R — siJ 2% \si R — siJ\s2 R — S2, 

si, S2 G C and R — Si, R — S2 G ^(3.8) 

The meaning of the integer R = s — r = ^P^ clearly results from 

[Na[P^IMI] = ^RMI (3.9) 
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We can now write the formula for the eigenvalues for the matrices A4{R): 
det {M{R)-fili) = 

R \ 
fi ■ 

- J 

g\R) 2si -R 2s2- R 



n 



det 



2tx 



Si{R - Si)jR - HSi{R - Si) \ S2{R - S2)\ R - IIS2{R - S2 



R 




9\R) 
27r 



E 



(2s - Rf 



s{R- s){R- fxs{R 



where we used the fact the a matrix || 6ij + ajOj || has all but one eigenvalues 
equal to 1 If we now exploit the symmetry s R — s of the C-invariant test 
vacua (where s E C <^==^ —seA has as the consequence that if s e C then 
—r = R — sEC)we get the eigenvalues 



R 



s{R 



and the equation (for R> 2 

271 



2 E 

seC,s<R~s 



seC,s<R — s 



2s -R 



(3.11) 



AR 



R^ 



s(R~s) 



/i 



(3.12) 



It is now immediate to see that none of the non trivial \A > can be the 
true vacuum since from ( |3.11| ) we find that there are tachyonic mesons with 



mass M^^^ 



-M 



2 (s-ry 



when r, s > or r, s < 0. 



Let us now consider the case of the trivial vacuum. In eq. ( |3.11| ) 

+ I and the eigenvectors 



the admissible values of s are now s 
associated to these eigenvalues are 



\^s{R) >-- 



6. 



s—r,R 



(Mt,-MV,,)|0> 



(3.13) 



and describe C-invariant states made of free particles. Similarly the sum- 
mation condition in eq. ( |3.12| ) becomes < s < -| . From the graph of 



(3.10) 
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the r.h.s. of ( |3.12| ) as a function of fi and from ^ (^^y) < it is clear 

that this equation always has ^ solutions and that it only has one negative 
solution when 

2n 2 ^ (2s -Rf 
0<-^<-^ E \(r. A (3-14) 

9 i^) ^ o<s<iR/2] - ^) 

The only way to avoid to have mesonic tachyons is to impose 

^>2E- + (E-- E -) + 2(^-2) Vi?>2 (3.15) 

2 

which implies in the limit 1 , i? < < A 



that essentially coincides with ( ^]2| ) a part from the constant —3: it is 
probable that some inconsistencies (tachyonic poles) show up in the range > 
> —3 when we consider the correlation functions of some (composite) 
operators. 

This relation will turn out to be fundamental in order to prove the Lorentz 
invariance of the mesonic spectrum in the thermodynamic limit, since the 
continuum version of eq. ( |3.11| ) is the consistency equation for the bound 
state wave function. 



4 The 't Hooft equation for the mesonic spec- 
trum. 

The 't Hooft equation is a light cone form for the Bethe-Salpeter equation, 
the aim of which is to find the spectrum of the bound states. This can be 
easily obtained from eq. ( ), ( |3.12| ) and ( ): the spectrum is given 
by 

^meson = 2yi? ^^R) = M^R^^{R) (4.1) 

where /i(-R) is any of the solutions of the eq. ( |3.12| ). But actually we are 
interested in the spectrum in the thermodynamic limit, i.e. when i? — > cxo 
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with ^ = const as L ^ oo. This means we are interested in the solutions of 
( p.l2| ) that scale as If we set = and s = xR, we can rewrite ( p.l2| ) 
in thermodynamic limit as 

J!^ = 2r^rf.ii^ (4 2) 

g2[R) X(l-X) l-ax{l-x) ^ ■ ' 

where both sides are divergent but their divergences cancel exactly, and what 
is left it is independent of R and yields a Lorentz invariant result for the 
meson mass. This is analogous to what happens in QCD2 when the 't Hooft 
IR cutoff A (|]lOl) is used. Moreover this sheds light on the nature of the bound 
state wave function renormalization as done in (|TI|): it is an expression of 
the fact that the IR divergences in the thermodynamic limit can be seen as 
UV divergences in the box. 

We can also get the 't Hooft equation in a more conventional way; we 
define the mesonic state as 



\ip,R>=J2Ms)Mlj,jO> (4.3) 

s 

and we get the 't Hooft equation requiring it to be the eigenstate of the mass 
equation = 2P^P~ , explicitly 

g\R) (\ 1 




27r \x l-xJJ^ \y l-y, 

The solution of this equation 



V9(0.4) 



1 — 2t 1 

^i^) = -7, 7— K{R) (4.5) 

x{l — x) 1 _ ^'-^meson 

x{l-x) M2 

does not seem to have a Lorentz covariant spectrum, but the consistency 
equation for the factor K(R) = ^-^^ j\_ dy Q — ^(y) is again ( ^]2| 
). 
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5 Conclusion. 



In this work we considered the Gross-Neveu model and we quantized it using 
the hght-cone approach in order to understand better how renormahzation 
and dynamical symmetry breaking come into the play. We discovered that 
we cannot describe the massless Gross-Neveu model and hence we cannot 
describe the "dynamical symmetry breaking". We found also that in the 
massive case the running of the coupling constant is only due to normal 
ordering effects but we have not an explanation for this phenomenon. 

A way to describe both the massless and the massive case is to consider 
the Yukawa lagrangian 

- ^ 1 1 

and let /i oo. 

Another interesting point would be to discuss the kinks in this formalism. 
All these things are left for future work. 



Acknowledgement. I thank P. Di Vecchia for useful discussions. 

6 Appendix 1. 

In view of eq. ( p.l6| ) we want to compute the following contributions 



T{x) = iip.^P+ + ^+i^^)ix) = NToix) + T^{x;xi,yi)NA[M^i)-Myi)] + -- 



xi,yi 



oo „ m 

m=0 i^i,yi)l<i<m j = l 



.11 



where T = 0(1) and we neglect the non leading contribution to the T. Since 
T.s are the terms with the highest power, we only have to consider in ( |2.14| 
) the contributions that come from the normal order of the pair ■ il^+{y) 

where the U{N) indices are contracted because the other are suppressed in 
j^. With this observation in mind we can set up a set of relations for the 
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coefficients in ([!])• The actual relations are 

ro(a;) = -g^ f K{x - y)To{y) -"^ f K{x-y) (.2) 



and 



Ti{x;xi,yi) = -g"^ K{x - y)Ti{y; Xi,yi) 

•'y 

-i^dz\xi-yi) 



5{xi - x) ( ^ - /To(yi) ) + 5{yi - x) ( ^ - g'^To{xi) 



and more generally (notice that it is not necessary to explicitly symmetrize 
w.r.t. the Xi and yi indices since the product N_A^\^XiLii^+{.Xi) ■ i^+iVi) ] is 
already symmetric) 

Tm{x;{xi,yi)i<i<rn) = - g"^ K (x - y)Tm{y; {xi, yi)i<i<m) 

■'V 

{yj,Xj)l<j<m-l)) (.4) 

where we defined 

K{x -y) = ^ idz\x - y){'ilJ+{x) ■ i/j+iy) - ^+{y) ■ ^+{x)) (.5) 
A simple computation yields 



■ y — x 



^Mx)-My) = j:^-^ (-6) 



and 



^(--y) = E^E^?^ (-7) 



It is now easy to get from ( [2| ) and ( [3| ) eq. ( |2.17| ) . 
The solution for eq. ( [4| ) reads 



l<i<m 



where g^y = g^{ 11 — g'^K)^y. This formula shows that all the terms can be 
expressed in function of the running coupling constant at different momenta. 



(.3) 
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